High temperature fractional quantum Hall states 
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We show that a suitable combination of geometric frustration, ferromagnetism and spin-orbit 
interactions can give rise to nearly flat bands with a large bandgap and non-zero Chern number. 
Partial filling of the flat band can give rise to fractional quantum Hall states at high temperatures 
(maybe even room temperature). While the identification of material candidates with suitable 
parameters remains open, our work indicates intriguing directions for exploration and synthesis. 



Landau symmetry breaking[l, 2] has been the stan- 
dard theoretical concept in the classification of phases 
and transitions between them. However, this the- 
ory turned out insufficient when the fractional quan- 
tum Hall (FQH) state [3, 4] was discovered. These 
highly-entangled quantum states are not distinguished 
by their symmetries; instead they are characterized by 
new topological quantum numbers such as robust ground 
state degeneracy [5, 6] and robust non-Abelian Berry's 
phases [7] of the ground states [8]. The new kind of 
order revealed in these topological quantum numbers 
is named topological order. [8, 9] Recently it was real- 
ized that topological order can be interpreted as pat- 
terns of long-range quantum entanglement [10-1 2]. This 
entanglement has important applications for topologi- 
cal quantum computation: the robust ground state de- 
generacy can be used as quantum memory [13]. Frac- 
tional defects from the entangled states which carry frac- 
tional charges [4] and fractional statistics [14-16] (or non- 
Abelian statistics [17, 18]) can be used to perform fault 
tolerant quantum computation[19, 20]. 

At present, the highly entangled gapped phases in FQH 
systems [3, 4] are only realized at very low temperatures. 
In this paper, we present a proposal to realize these states 
at high temperatures (even room temperature). The 
ideal is to combine spin-orbit coupling, ferromagnetism, 
and geometric frustration. Both spin-orbit coupling and 
ferromagnetism can have high energy scales and can ap- 
pear at room temperature. In some cases, combining 
them leads to energy bands with non-zero Chern num- 
bers and filling such an energy band will give rise to 
integer quantum Hall states. Further, in geometrically 
frustrated systems — lattices on which hopping is frus- 
trated — some of these topologically non-trivial energy 
bands can be very flat [21, 22]. These would mimic Lan- 
dau levels in free space. When such a flat band with a 
non-zero Chern number is partially filled (such as 1/3 or 
1/2 filled), FQH states can appear. Here we study a sim- 
ple example of this idea on the geometrically frustrated 
kagome lattice. 

Several aspects of the above ideas have been studied 
actively in recent research. Spin-orbit coupling can lead 
to a topological insulator in various geometrically frus- 




FIG. 1: The kagome lattice is a triangular Bravais lattice 
with a 3-point basis labelled / = 1,2,3; ai = x and a2 = 
(x + a/3^)/2 are the basis vectors. In the metallic kagome 
lattice Fe3Sn2, spin-orbit coupling arises from the electric 
field due to the Sn ion at the center of the hexagon. 

trated systems, [23-26] and non-collinear magnetic order 
can lead to integer quantum Hall states. [24, 27] Alter- 
natively, interactions in geometrically frustrated systems 
can break time-reversal symmetry [28-36] which again 
can give rise to integer quantum Hall states. Here we 
will show that a natural extension of these ideas may 
set the stage for FQH states and other highly entangled 
states with fractional statistics and fractional charges — 
possibly even at room temperature. 

Nearly flat band with non-zero Chern number We con- 
sider nearest and next-nearest neighbor hopping on a 
kagome lattice with spin-orbit interactions described by 
the following Hamiltonian 

{ij)a 

where creates an electron with spin a on the site 
on the kagome lattice. Here {ij) denotes nearest neigh- 
bors and {{ij)) next-nearest neighbors. The second and 
fourth terms describe spin-orbit interactions which pre- 
serve time-reversal invariance. R^j is the distance vector 
between sites i and j and Eij the electric field from neigh- 
boring ions experienced along R^j. 

To obtain FQH states, we need to break time- 
reversal symmetry. This is likely to happen sponta- 
neously from exchange effects in the flat band that cause 
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ferromagnetism [28-30]. Alternatively, one can apply an 
external magnetic field or couple the system to a ferro- 
magnet. In the extreme limit the electron spins are to- 
tally polarized within the partially filled band, which is 
the case we examine here. Hence we consider spin-orbit 
coupling that also conserves Sz^ i.e. the electric field on 
each site is in the 2D plane. 

We first study just nearest-neighbor hopping so t2 = 
A2 = 0. In momentum space, Eq. 1 becomes 



i^k = -2ti 
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where ai = f , a2 = + a/3^)/2, a3 = a2 — ai and 
= k • a^. We use units where the hopping parameter 
ti = 1. The +(— ) sign refers to spin up (down) electrons; 
from here we focus on just the spin up electrons. 

The spectrum consists of three energy bands and is 
gapless at Ai = 0, ±a/3. At all other points the spectrum 
is gapped and the top and bottom bands have unit Chern 
number with opposite sign while the middle band has 
zero Chern number. Here the Chern number is defined 
as [37] 



^ / d'kFi2{k) 

^TT Jbz 
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where Fi2{k) is the associated field strength given by 
^12 (^) = ^A2(k) - ^Ai{k) with the Berry connec- 
tion Aij^{k) = — ^(nk|^f-|nk). In the above |nk) is a nor- 
malized wave function of the respective band. 

Focusing on the lowest band which has a non-zero 
Chern number, we look for a regime in which this band is 
very fiat compared to the bandgap and the energy scale 
of interactions. We denote W as the maximum band- 
width of the lowest band, A12 as the minimum bandgap 
between the two lowest bands and U as the strength 
of electron-electron interactions in our system. When 
U interaction effects dominate kinetic energy and 

partially filling the fiat band would favor the Laughlin 
state. [4] Since band mixing could destroy the band fiat- 
ness, ideally A 12 ^ U. Hence we aim to maximize the 
ratio Ai2/I^ in order to obtain FQH states. 

As the middle band has zero Chern number, any mix- 
ing between only two of the bands would not change the 
Chern number of the lowest band. If the lowest band 
remains flat even with mixing then A13, the minimum 
bandgap between the lowest and highest bands, (and con- 
sequently the ratio Ais/W) is also of interest. 

We find that W > A12 always — as the bandwidth 
vanishes, the bandgap between the two lower bands van- 
ishes in the same way (due to topological symmetry in 




FIG. 2: (Color online) Results for nearest neighbor hopping 
as a function of Ai (nearest neighbor spin-orbit coupling). 
The bandwidth of the lowest band W vanishes at Ai = a/3; 
however W > A12 always where A 12 is the bandgap between 
the two lowest bands. Here we show the band structure for 
Ai = 1 where W 1.3 and A13/W 3.7. The spectrum 
does not have a clear separation of energy scales between A 12, 
W and U, the interaction strength, which could be due to 
limitations of the nearest neighbor hopping model. 



real-space) [38] , see Fig. 2. Here we show the band struc- 
ture for Ai = 1 where W 1.3 and Ai3/VI^ 3.7. As 
A12/W < 1 always, the spectrum does not have a clear 
separation of energy scales. When interactions are on the 
order of the bands will mix. This scenario is quite 
different from Landau levels in free space that are fiat 
and well-separated, which could be due to limitations of 
this simplest model. 

For a more realistic scenario, next we include second- 
nearest neighbor hopping. This gives us additional terms 
in the Hamiltonian 
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In this larger parameter space the band maxima and 
minima are no longer fixed at the same symmetry points. 
We find that the largest values of Ai2/1^ (and A13/VF) 
occur when Ai and A2 are of the same sign - in which 
case the results are symmetric under changing signs of 
both A's in this spin polarized case. In Fig. 3, A12/W 
is plotted as a function of t2 for three values of Ai = 
A2. We see that at negative values of t2 a lower spin- 
orbit coupling is needed; while for positive values of t2 
higher values of spin-orbit coupling would maximize the 
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FIG. 3: (Color online) With the inclusion of next-nearest 
neighbor hopping, we obtain much higher bandgap to band- 
width ratios. We choose three values of Ai = A2 = 0.3, 0.5 
and 0.7, and sweep A12/VK with t2- For lower values of spin- 
orbit coupling, the ratio peaks at negative t2; for relatively 
higher values of spin-orbit coupling the converse is true. 




FIG. 5: (Color online) Three fairly flat bands are mutually 
well-separated with bandgap to bandwidth ratios of A12/VK = 
8.7 and A12/VK = 24. Parameters used are t2 = — 0.3, Ai = 
0.6 and A2 = (Case 2). 




FIG. 6: (Color online) Distribution of the field strength Fi2{k) 
(Eq. 3) in the Brillouin zone for the flat bands in Cases 1 and 
2 discussed above. They do not contain sharp features — 
especially Case 1 with the flatter band — hence the magnetic 
length scale remains on the order of the lattice constant a. 



FIG. 4: (Color online) A very flat lowest band — well- 
separated from the two higher bands — is obtained with the 
parameters t2 = —0.3, Ai = 0.28 and A2 = 0.2 (Case 1). 
The bandgap to bandwidth ratios are high: A12/W = 52 and 
A13/W = 99 respectively. 



bandgap to bandwidth ratio. 

We present two examples where A12/W and Ais/W 
reach high values at ^2 = — 0.3. In Case 1, setting ^2 = 
—0.3, Ai = 0.28 and A2 = 0.2, we obtain a very flat lowest 
band separated from the two higher bands by a large gap 
(see Fig. 4). The values of Ai2/I^ and A13/W are 52 
and 99 respectively. In another example, all three bands 
are fairly flat (particularly the lowest one) and mutually 
well-separated. The parameters used are ^2 = —0.3, Ai = 
0.6 and A2 = 0. In Case 2, we obtain A12/W = 8.7 and 
Ais/W = 24 respectively (see Fig. 5). 

Calculating the Chern number c of the lowest flat band 
in these two cases, we find it is 1. This is expected as 
slowly turning off ^2 and A2 does not close the bandgap 
— and we have previously seen that in the absence of 
next nearest-neighbor hopping, the lowest band always 
has unit Chern number. When A ^ U ^ W is satisfied, 
partial filling of this flat band would favor the FQH state. 

The distribution of the field strength Fi2{k) in the Bril- 



louin zone is plotted in Fig. 6. We observe there are 
no singularities or very sharp features but Fi2{k) varies 
fairly smoothly especially in the first case with the flat- 
ter band. The presence of singularities — e.g. localized 
at the Dirac point — would have signalled a new (and 
much larger) length scale in the system. In our case, both 
the magnetic length scale (arising from spin-orbit inter- 
actions) and the variation of field strength Fi2{k) are on 
the order of the lattice constant a. 

Thus the interaction energy scale is generated from 
the lattice constant a: U ^ e^/ea where e is the di- 
electric constant. If the flat band is 1/3 filled, analo- 
gous to results from FQH states in semiconductor-based 
systems [39] the gap for this v = !/?> state is here roughly 
0.09e^/ea = 500K. (We choose e = 3 and a = lOA; a is 
defined as the square root of the unit cell area). As the 
interaction energy scale is a hundred times larger than 
in semiconductors, we may see the FQH effect at room 
temperature. As band gaps A12 and A13 are easily much 
higher than room temperature, a fully filled band could 
give the integer quantum Hall state at unusually high 
temperatures too. 

Materials realization We see that a suitable combina- 
tion of geometric frustration, ferromagnetism and spin- 
orbit interactions can give rise to nearly flat bands with 
a large bandgap and non-zero Chern number. The fer- 
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romagnetism can arise from an external magnetic field, 
a ferromagnetic substrate for a thin film sample, and/or 
exchange effects. If the flat band is close to the Fermi en- 
ergy, partial filling of the flat band can be controlled by 
doping and give rise to FQH states at high temperatures. 

The choice of parameters in our calculations is based on 
known values of spin-orbit coupling. For example, in Her- 
bertsmithite (a common copper-based 2D kagome lattice) 
the spin-orbit interaction is 8% of the kinetic energy. [40] 
Other compounds with 4d or 5d orbitals (instead of 3d as 
in Cu) may experience a larger spin-orbit interaction. For 
instance, the strength of spin-orbit coupling in Iridium- 
based kagome compounds can be on the order of mag- 
nitude of the kinetic energy, hence the substitution with 
4d or 5d atoms in metallic kagome lattices could result 
in hopping parameters similar to the ones used in our 
work. Alternatively, making thin films of frustrated lat- 
tices with 4d or 5d atoms may lead to a flat band with 
strong spin-orbit coupling, where exchange effects in this 
flat band could cause ferromagnetism. 

Most existing kagome compounds are Cu-based insu- 
lators. Some 2D kagome lattices show metallic behavior, 
for instance Fe35'n2.[41, 42] This material shows ferro- 
magnetism along the c-axis above 60K and in the kagome 
plane below this temperature. Also, as spin-orbit inter- 
actions can be simulated in cold atom systems, [43, 44] it 
is possible to realize our hopping model in such systems 
as well. This would provide a method to obtain FQH 
states in cold atom systems. 

In short, flat bands with non-zero Chern number arise 
in the examples we have given and in other geometri- 
cally frustrated systems with suitable levels of spin orbit 
interaction. By partially filling these bands e.g. via dop- 
ing, one can expect the emergence of FQH states at high 
temperatures. While the identification of exact material 
candidates remains open, our work indicates intriguing 
directions for synthesis and development. 

After completion of this paper, we learned that T. Ne- 
upert et. al[45] also discussed the possibility of the FQH 
effect in interacting two-band lattice systems, while K. 
Sun et. al[46] found flat bands with non-zero Chern num- 
bers on various lattices (including the kagome lattice) 
after including some complex hopping. 

This research is supported by NSF Grant No. DMR- 
1005541. 
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